DARBOUX TRANSFORMATION FOR 
CLASSICAL ACOUSTIC SPECTRAL PROBLEM 
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Cl^l Abstract. We study discrete isospectral symmetries for the classical acoustic spectral 

' problem in spatial dimensions one and two, by developing a Darboux (Moutard) trans- 

formation formalism for this problem. The procedure follows the steps, similar to those 
CO ' for the Schrodinger operator. However, there is no one-to-one correspondence between the 

two problems. The technique developed enables one to construct new families of integrable 
potentials for the acoustic problem, in addition to those already known. 

The acoustic problem produces a non-linear Harry Dym PDE. Using the technique, we 
reproduce a pair of simple soliton solutions of this equation. These solutions are further 
Ph I used to construct a new positon solution for this PDE. Furthermore, using the dressing 

chain approach, we build a modified Harry Dym equation together with its LA-pair. 

As an application, we construct some singular and non-singular integrable potentials 
(dielectric permitivity) for the Maxwell equations in a 2D inhomogeneous medium. 



Introduction 



> 

o 

^ I This note develops the Darboiix transformation and dressing chain formahsm for the 

[ classical acoustic spectral problem (below just the "acoustic problem") and the related 
Harry Dym (HD) equation. It treats the problem in the same vein as it is done for the 
Schrodinger operator and the related KdV (mKdV) hierarchies. The acoustic problem 
I and the Schrodinger operator are closely connected. This connection constituted a base 
[ for the approach to the acoustic problem and the HD equation in [1-3]. However, as is 
discussed below, the relation between the problems is not utterly straightforward. 

The acoustic problem describes wave propagation in inhomogeneous acoustic or elec- 
tromagnetic media and just like the Schrodinger equation is non-integrable for an ar- 
bitrary potential. For applications, it is important to be able to construct intergable 
potentials, which result in solutions with given properties or asymptotic behavior. For 
instance, for the purposes of transmission of information, reflexionless potentials are im- 
portant. These potentials are such that the problem admits solutions, which asymptote 
to e~^P^ as X — oo and T{p)e^P^ as x — >• oo, with the passage coefficient T{p) G C 
being one in absolute value. This was recently studied in a work by Novikov [4], which 
has drawn our attention to the problem. 

The latter work constructs a family of so-called B-potentials for the acoustic problem 
via a semi-classical solution ansatz. We show that these potentials naturally come up as 
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a result of a one-step Darboux transform on the "vacuum background". In continuous 
media electrodynamics, potentials of the acoustic problem can be interpreted as the 
medium's dielectric permitivity. The B-potentials in question possess point singularities, 
and hence their physical meaning is not entirely clear. On the other hand, the dressing 
technique enables one to construct a variety of potentials, which are neither B-potentials, 
nor singular. We illustrate it by a single act of dressing chain closing (in dimension one 
and two) which yields regular integrable potentials (dielectric permitivity). 
The acoustic problem on the real line is described by the following ODE: 

A 

^xx = -^T^^- (1) 

This equation models wave propagation in non-homogeneous (acoustic or electromag- 
netic) media. Consider for instance the Maxwell equations in a medium without exter- 
nal sources with the standard notations (E, H) for the electromagnetic field, as well as 
D = eE, B = kH. Suppose, the medium is isotropic but inhomogeneous with k = 1 and 
e = e{x, y, z). Then 

1 c)D 1 9B 
rotB=- — , rotE = — , divD = divB = 0. (2) 

c ot c ot 

One can easily exclude the quantity B from (2) to obtain an equation connecting the 
quantities E and D: 

1 a^D 

"t"tE=-jjg^. (3) 

For the electric field E one writes E = e*'*'*i/?, with i/? = tl){x^ y, z) and taking into account 
the last equation of (2) obtains: 

V((3^)+A^ = -^.^, (4) 

where A is a three dimensional Laplassian. 

The equation (1) follows if one lets e = e(a;), = (0, 0, ip{x)), A = —uP'jf?^ u~'^(x) = 
e{x). The dielectric permitivity e{x) will be henceforth referred to as a potential. 

Alternatively, one can choose e = €{x,y) as well as if) = {0,0,i(^{x,y)). If this is the 
case (4) is reduced to a linear PDE 

All; = Xeip, (5) 

where A is a two dimensional Laplassian. 

Hence, studying the equations (1,5) is of interest for continuous media electrodynam- 
ics. A similar case can be made in acoustics, whence comes the original name of the 
equation (1). In both cases the equations describe transmission of signals, and are quite 
relevant for applications. 

For the Schrodinger equation one of the most efficient ways of building potentials 
allowing exact solutions, is the method of factorization, or the Darboux transformation. 
Developing a similar formalism for the equation (1) appears to be a natural thing to do. 
Its basics are presented in the next two sections of this article. In particular, we derive 
the related generalized Crum formulae, build chains of discrete symmetries and study 
their simple closing. We argue that there is no one-to-one correspondence between this 
formalism and the well-known technique for the Schrodinger equation [5,6,10]. 
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The acoustic problem (1) is also interesting from the integrable systems viewpoint. It 
is known that this equation represents the L-equation of the Lax pair (or the LA-pair) 
for the nonlinear Harry Dym (HD) PDE: see [1-3] and references therein. We illustrate 
the Darboux transformation technique for the ODE (1) by constructing a new positon 
solution of the HD equation from a pair of its simple soliton solutions. 

The Darboux transformation is known not only as a way of finding exact solutions 
of nonlinear equations, but also as a resource for proliferation of these equations, that 
is building new integrable PDEs together with their Lax pairs by means of the dressing 
chain technique. As an illustration, we construct a modified Harry Dym (mHD) equation 
(and its Lax pair) which has a remarkably simple form. 

In the last section we turn to the equation (5). The Moutard transformations for 
this equation provide a simple method for construction of exact solutions of the Maxwell 
equations (2) with a dielectric permitivity e = e(x, y), which in general has a complicated 
singularity structure. However, a simple periodic closing of a dressing chain, generated 
by the Moutard transforms results in a regular integrable 2D dielectric permitivity. 

Discrete symmetries of the one-dimensional acoustic problem 

The Darboux transformation technique makes use of the existence of specific discrete 
isospectral symmetries of the equation under consideration and is standard in the theory 
of integrable PDEs: e.g. [5,6] and references therein. 

We start out with a standard non-linear substitution [1,2,4] due to which typically the 
solutions of (1) will be represented parametrically: 

u{x)=Vy{y), x = v{y). (6) 

This reduces the equation (1) to 

ll;yy = UlPy+Xlp, (7) 

with U = Vyy/vy. Thc above quantity U{y) will be referred to as a potential as well as 
the dielectric permitivity function e{x) = ^^^^ mentioned earlier. 

In spite of the fact that (7) is easily reducible to the Shrodinger operator, from the 
point of view of finding integrable potentials this connection is not trivial. Let us address 
this issue in more detail. 

A substitution 

transforms (7) into the stationary Schrodinger equation 

^yy={X + V{y))i;, (7.1) 
where the potential V{y) is related to the potential U (y) of the acoustic problem (7) via 

-2Uy , , 

V= ^. 7.2 

4 

Linearizing (7.2) with a substitution U = —2py/p, one sees that p{y) in turn satisfies 
(7.1) with A = 0. 

The Darboux transformation for (7.1) is well known and at the first sight it may appear 
that developing an independent technique for the acoustic problem (7) is superfluous. 
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However, the following argument shows that this is not the case. Namely, there is no 
one-to-one correspondence between the problems (7) and (7.1). 

Indeed, let U{y) be a specific potential for the acoustic problem (7), not depending 
on any free parameters. From (7.2) one can get (uniquely) the Shrodinger potential 
V{y) and further substitute it into (7.1). In order to reconstruct the initial potential 
[/(y), equation (7.1) should be solved with A = 0. Let the solution he p = p{y, Ci, C2), 
depending on a pair of constants Ci,C2. One of these constants, say Ci, plays the 
normalizing role and can be omitted. Then the restored potential U = U{y, C2) will 
depend not only on 2/, but the free parameter C2 as well. Hence, a single potential 
in the Shrodinger operator generates the whole family of potentials for the acoustic 
problem, and in order to single out a specific potential for the latter one would have 
to subsequently develop some selection mechanism by studying the sequence of maps 
Uiij) V{y) t/(y, C2) U{ii). This necessity gets bypassed if one develops the 
Darboux transform formalism directly apropos of the operator (7) without using (7.1), 
and this is done in this and the following section. 

Following Shabat [5] we seek elementary discrete symmetries of the equation (7) ef- 
fecting the change 

^^^(1)=/^^ + ^^, (8) 

for some A-independent functions / and g of y. 

One easily verifies that there are three distinct discrete symmetries of the type (8) for 
the equation (7). They are 

^^^i') = t, v^vi^) = L ^^^w = ty^ y^yW=f^ (9) 

V V Vy J Vy 

and 



Vy ^ vi'^ = Vy {^) , U^U(^^ = U + 2Dln^ 



(10) 



In the latter equation '0i = 'il;i{y,Xi) is a particular solution of (7) with the spectral 
parameter value Ai, further referred to as a "prop solution", D = dy and i/^i^y = -DV'i- 

Note that the former two symmetries (9) define the new quantity ' ~ -j^ in a 

way independent of any solution '0(|/, A) of (7). These symmetries arise as a particular 
case of (8) as the result of gauging corresponding to the choice oi f ov g alternatively 
zero. These symmetries have a trivial kernel in the solution space of (7). According to 
the terminology of [5] we call the symmetries (9) T-symmetries, sometimes referred to 
as Schlesinger transforms^. 

On the other hand, the transformation (10) alias the Darboux transformation, which 
[5] calls an S-symmetry, does have a non-trivial kernel on the solution space of (7) (one 
can let ip = ipi in the first equation of (10) and get zero). This property will be essential 
in the sequel. 

Along the way, we shall use the popular term "dressing" for the application procedure 
of the transformation (10) to a triple (ip^v^U), the resulting pair (il;^^\v^^\U^^^) being 
referred to as the "dressed" one. 



^In the context of soliton solutions, the T-symmetries play the part of explicitly invertible Backlund 
transforms [7]. 
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Despite a nearly trivial countenance, the Darboux transform (10) has a remarkable 
capacity to enable one to engineer potentials with arbitrary discrete spectra ad hoc. 
Indeed, suppose it is possible to solve the equation (7) formally (namely, obtaining among 
others some "non-physical" solutions which are not in L^) for some potential U and all 
A e IR. Suppose, 'ipi{y,Xi) is such a solution. Let us denote its linearly independent 
counterpart as i/jid/, Xi), i.e. 



Dressing ipi according to (10), we find 

7(1) 

Therefore, if one comes up with a non-physical prop solution t(ji by requiring that its 
derivative ifji^y be strictly positive and rapidly growing as \y\ — > oo, then in the spectrum 
of the dressed potential U^^^ there will appear a level Ai, not present in the original 
spectrum for U. Since the principle for the choice of the value of Ai is such that this 
value is not to be present in the physical spectrum for U, repeating the dressing procedure 
n times will result in a potential t/^"^ possessing n new pre-chosen levels Xj, j = 1, . . . ,n. 
See (13,14) in the sequel. 

Conversely, the function generates an inverse transformation (undressing) to (10). 
Thus, one can as well remove some pre-chosen levels from the spectrum of a potential. 

Crum formulae and dressing chains for 

THE classical ACOUSTIC SPECTRAL PROBLEM 

Below we present the formulae describing an n-step dressing procedure for any n > 1, 
whose analogues are known for the Schrodinger equation as the Crum formulae [8]. We 
derive them for equation (7) following the procedure exposed in [9]. 

A single act of dressing (10) can be iterated n times to yield a triple v^'^\ C/*^"^). 

One starts out by dressing a triple (V', f , U) = v^^\ U^^^) corresponding to a spectral 

parameter A with a prop function V'l = V'l^^ which is a formal solution of (7) with a 
spectral parameter Ai and the potential U^^\ The resulting solution i/j'^^'' solves (7) 
with the dressed potential U^^^ (and the same spectral parameter A). On the jth step, 
j = 1, ... ,n one uses some prop solution ipj'' which solves (7) with a pre-dressed 
potential U^^~^'^ and a spectral parameter value Xj to produce the j times dressed solution 
t/j'^^^ and potential U^^'' (as well as the function v'^^^ with U^^'' = VyJ /vy '^). Note that the 
spectral parameter A in the dressed equations for ip^^^ is the same for all j = 1, . . . , n. 

It's easy to see that the n times dressed solution V'^"^ shall have the form 

n 

= J2 djDH + (-1)>, (11) 

with the functions-coefficients ttj to be found, which of course will depend on the choice 
of the prop solutions ip^^ ^\ It follows from (10) that 



= U + 2D\Yia, 
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(12) 



for 



where the eUipses in the first formula stand for the terms containing the derivatives of 

if)'^^^ of orders from 1 through n — 1. 



So far the choice of the prop solutions i/)^"^ has been quite arbitrary. But suppose 
now that the original equation (7) possesses n distinct formal solutions corresponding 

to spectral parameter values Xj, j — 1, . . . ,n. Let ipj = ip^^'^ and consider the following 
dressing procedure (which will be further used for the dressing chain construction): 



^(1) 

^(n-l) 



(0) 












(0) 








€1. 








(0) 







t/(0) 
[/(n-i) 



(*) 



Namely, for j = 1, . . . , n on the above diagram (*) every new line j + 1 is obtained by 
dressing the functions from the preceding line j by (10) with a prop solution i/jj'' """^ 
marked in bold. 

Zeroes, proliferating as one moves down the diagram stem from the non-trivial kernel 
property of the S-symmetry, and it is this property that now enables one to find the 
unknown functions aj. Indeed, substitution of any il^j = ip^^^ for ip in the right hand 
side of (11) shall yield zero. Hence, the coefficients aj satisfy a system of n independent 
linear algebraic equations, namely 



^afeDV, + (-l)>, =0, j = l. 



n. 



fc=i 



Solving it by the Kramer rule and substituting the result into (11) and (12), we end up 
having: 



= t/ + 2L»ln-^, 



(13) 



where A„, A„ are determinants of square nxn matrices, whereas A„_|_i - of an n+l xn+l 
matrix as follows: 



An = 



An = 



(14) 



^n - D^ll^n 
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Note that the hnearity of (7) makes the choice of the sign before ip'^'^^ irrelevant. 

The obtained formulae (13,14) make it possible to find rich families of exact solutions 
of (1). The easiest case is dressing from the birthday suit, or on the vacuum background, 
assuming U — (hence x — ciy + C2, u{x) — ci, where ci and C2 are arbitrary constants). 

Such a natural rigging yields gratis all the B-potentials reported in [4] , moreover the 
formulae for their computation derived therein turn out to be particular cases of (13,14) 
with merely Vy = 1. For instance tfji = sinh^i, il)2 = cosh ^2 (here n = 2, yj are constants, 

= kj{y — yj), j — 1,2) yield a reflexionless B-potential with a power 2/3 singularity: 
f(^) = — — wherp^ 



\ 2 

/c2 sinh ^1 sinh ^2 — ki cosh ^1 cosh ^2 



ki sinh ^1 sinh ^2 — ^2 cosh ^1 cosh ^2 
The same i/ji and '02 = sinh ^2 yield another potential with a power 4/5 singularity: 
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(-2) / /c2 sinh^i cosh ^2 — ki sinh ^2 cosh^i 
\ ki sinh ^1 cosh ^2 — ^2 sinh ^2 cosh ^1 

By construction, these potentials have only two levels Ai^2- 

In addition, all the regular reflexionless potentials can be also built by the formulae 
(13,14) once again by dressing U = 0. Here comes the proof. The passage coefficient for 
a regular n-level reflexionless potential can be expressed by a well-known formula 



Tn{p) = n 



kj — ip 

\A ' 

with Xj — kj. As has been pointed out earlier, the levels Xj can be successively removed 
from the spectrum by means of the inverse of the Darboux transformation (still having 
the form (10)), each application of which will kill a term in the product. Successively 
applying this procedure n times, for the passage coefficients we have 

Tn{p) Tn-i{p) ... To{p) = 1. 

This proves our assertion, because the case R{p) = and T{p) — 1 for the reflection and 
the passage coefficient for all p is feasible with U = only. It's worth reiterating the 
point that the argument above owes itself to the fact that the S-symmetry possesses a 
non-trivial kernel in the space of solutions of (7). 

Matveev and Salle [6] flnd super-reflexionless potentials for the KdV equation, alias 
positons. In the same vein one can operate on the equations (1,7). In order to do 
so, one should use the formulae (13,14) with n = 2 choosing the prop solutions 1(^1^2 
respectively as 'i/'i(y, Ai) and 'i/'id/, Ai + 5), and then letting 5 ^ 0. If withal U = and i/^i 
generates a single soliton potential, then (13) deflnes a single positon potential. The next 
section describing the positon solutions of the HD equation contains the aforementioned 
computation. 



^Note that the expressions for m^^' are parametric. In order to interpret the formulae correctly, the 
reader is referred back to (1,6,7). The orders of the singularities pertain to the potential e{x), which is 
a zero of the function u{x) and a singularity of the potential U{y) — Vyy/vy, where the function v{y) 
solves the equation u[v{y)] = Vy{y). 
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In addition to B-potentials, various other interesting ones can be produced. For in- 
stance, one can construct soluble potentials with a finite equidistant spectrum^. 

One can also investigate potentials which change in a specific simple way under the 
Darboux transform, e.g. such that U — > CZ+const or U — >constC/. For the Schrodinger 
equation, the former transformation is shape-invariant and for n = 1 results in the 
harmonic oscillator potential. Let's develop an analogue for the model (1,7) under inves- 
tigation. Suppose, 

for a constant u. Then we can obtain parametrically the function u{x) from (1) as 
follows: 

u{x) = 0:2; exp ^—(jj^z ^ , X = xq — auP' J dzexp ^—uj'^z ^ , 

where, k, Xq, ol arc real constants and z = exp(— y/cu^). The prop function ifji rendering 
the potential U^^'> from U has the countenance = exp(— cu^z) and solves (7) with an 
eigenvalue Ai = It's easy to verify that the dielectric permitivity e{x) — l/u^{x) 

has a second order pole at x = xq. 

The theory of the Darboux transformation for the Schrodinger equation utilizes the 
concept of dressing chains of discrete symmetries and their closing. The work of Veselov 
and Shabat [10] elucidates how the dressing chain closing method can be used in order 
to obtain various potentials with meaningful mathematical physics. Namely, a simple 
closing procedure leads one to the harmonic oscillator potential (resulting also in a shape- 
invariant change of potential). A more complicated closing results in finite-gap potentials 
as well as the fourth and the fifth Peinleve equations, see [10]. 

Dressing chains can be written out for the equation (7) as well. Let us introduce a 
sequence {/n}n>i of functions as follows: 

with the quantity i/'rT as it has been introduced in the diagram (*) above (where 
it appeared in bold). In particular, it corresponds to the pre-chosen value An of the 
spectral parameter. 

One can verify by hand starting from n = 1 that 

n 

Besides, direct substitution shows that fn satisfies the equation 

fn + fn-U^''^fn = K, 

where /' = Df. The two latter relations imply the recursion connecting and fn+i as 
follows: 

ifnfn+lY = fnfn+lifn ~ fn+l) + ^n+lfn ~ ^nfn+1- (15) 

This equation (15) represents a dressing chain for the acoustic problem. 



^The same statement applies to the (stationary) Schrodinger equation. 
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In a way analogous to the theory of dressing chains for the Schrodinger equation [10] 
wc arc interested in T-periodic chain closing, namely imposing the condition fn+T = fn 
for an integer T > 1. We shall consider here the easiest case T = 1. 

Given the spectral parameter values Ai,2, one obtains a one-parameter family of po- 
tentials, indexed by a constant c: 

^ -(^1 - ^2)'y' + 2c(A2 - Xi)y + 6Ai - 2A2 - 
2[(Ai-A2)y + c] 

If A2 = 3Ai > and c = 0, we can express the function u{x) parametrically: 



x{y) = — Erf(a|/), u{y) = exp {-a'^y) , 

where = -A/2 > 0. 

It is known that for the Schrodinger equation, a nontrivial chain closing operation 
with T > 1 results in finite gap potentials [10]. Such potentials for the HD equation are 
due to Dmitrieva [3] . A close analogy between the Schrodinger equation and the acoustic 
problem [1,2,11] suggests that one can expect results similar to those of [10] apropos of 
the analysis of higher order chain closing for T > 1. We expect that potentials built in 
such a way can have interesting physical applications, such as for instance a model of 
wave propagation in media whose dielectric permitivity is a periodic function of a single 
spatial variable. 



HD AND mHD equations 

The 1-^1 HD equation 

Ut = U^Uxxx + Pux, (16) 

with some real constant (3, has been studied quite extensively since late 70's: see [1-3] 
and references therein. It arises in the study of evolution equations solvable via the 
spectral transforms method based on the string rather than the Schrodinger equation. 
The principal approach to it has been based on its relation to the KdV, mKdV and other 
more classical hierarchies of integrable PDEs [2,3]. However, as was shown above, this 
relation is not entirely straightforward, and the direct approach developed herein enables 
one to produce new solutions of the HD equation in addition to those already known. 
As an example, below we construct a simple positon solution. 

The acoustic problem (1) is the first equation in the LA-pair for the HD equation (16), 
the full pair being 

(17) 

i^t = (4Au + /3) V'x - 2Xu^'il;. 
The coordinate change (6) in the presence of time dependence becomes 



thus 



dx —dy, dt^dt dy. 

Vy Vy 
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After this change (17) becomes 



y 



(18) 

4a) ^l:y - ^V', 



and the HD equation (16) transforms to 

Vy{VtVyy - VytVy) + 3V^j^ + Vy{V4.yVy - Av^yVyy + (3Vyy) = 0. (l9) 

with the notations v^y, v^^y for the partial derivatives in y of order 3 and 4 respectively. 

The goal now is to extend the Darboux transformation (10) for the equation (7) alias 
the first equation in (18), so that it agrees with the second equation in the Lax pair. One 
just includes the t-dependencies in (10). At this point it only provides the value of the 
partial derivative 
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rather than the dressed quantity v^^\y^t) of interest. Hence, let v^^^ = B{y^t) be un- 
known and let's assume that v'^^^ satisfies the second equation of the pair (18) with the 
dressed according to (10) function the quantities (A,/3) remaining the same. One 
can express the unknown quantity B as follows: 

B=(^\ ((3 + 4i,Vy + vt-2vyUy) + ^^l^-4vy-(3, 

\'Vl,y / Vl,y 
and verify that By = At. It follows that 

v^^\y,t) = j Ady + Bdt, (20) 

with a closed 1-form under the integral. 

Hence, the Darboux transformation (10,20) is an S-symmetry for the LA-pair (18), 
and therefore for the HD equation (19). 

This enables one to construct exact solutions for this equation. We exemplify it with 
a single soliton solution and a single positon solution. 

Let V = y, Xi = i(^i{y,t, k) = sinh[0(?/, t, k)] with (f) = k (y + (4A;^ + P)t), then by 
(20) one has 

= J_(0_tanh</.)-(4 + /3)t. 

The function v^^\y,t^k) determines a single soliton B-potential U^^'> = Vy]j /vy^\ men- 
tioned in the previous section. 

The single positon potential is obtained from two distinct soliton solutions ilji{y,t, k) 
and ilji{y,t, k + S), using them as the prop functions ■01,2 in the formulae (13,14) with 
n = 2 and taking the limit as 5 — > 0. Namely, 



"y "y 



''Pl,yyk'^l,y ~ '^l,yy'^l,yk 
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where the subscript ^ means differentiation by k. Taking ij^i exphcitly as the hyperbohc 
sine in the previous example results in 

^(2) / sinh(20) + 20 Y 
^ \sinh(20) - 20y ' 

with (f)^k{y + {I2k'^ + I3)t). 

It is well known that the dressing formalism enables one to produce hierarchies of inte- 
grable PDEs. Borisov and Zykov [12] proposed a technique for proliferation of integrable 
equations, which they applied to the KdV and the Sine-Gordon (SG) equations. The 
technique is based on the discrete symmetries dressing chain closing. The main idea of 
the approach is as follows. The equation (for illustration purposes let us take the KdV 
equation) is written as a compatibility condition of a pair of equations, further denoted as 
Li and A^. Each of these equations is quadratic in the auxiliary field. Using invariance 
of the pair with respect to the Darboux transformation (which is viewed as a discrete 
symmetry), a second pair L2, A2 of equations is built. Excluding the potentials from 
Li, L2 and Ai, A2, it is possible to obtain two equations, which [12] calls an x and a t 
chain, respectively. (We further use the notations Cx and Ct instead.) If a potential is 
excluded from Li and ^1, one ends up with a modified equation mKdV. The equations 
Cx and Ct can be converted into the Lax pair for the mKdV equation in two ways, the 
Darboux transformation being already known. 

This procedure can be repeated, producing new equations with their LA-pairs. In 
this vein, the equations m^KdV and m^KdV were obtained. The former becomes the 
exponential Calogero-Degasperis equation [11] after an exponential change, the latter 
contains an elliptic equation of the same authors. 

In spite of its simplicity, the technique described is very powerful. This can be illus- 
trated by the following examples. First, see [13], the m-^KdV equations with = 0, 3 
together with the Krichever-Novikov equation exhaust (modulo a contact transforma- 
tion) all the integrable equations of the form ut + Uxxx + f {uxx^Uxiu) = 0. Second, 
applying their approach to the SG equation, the authors of [12] have succeeded to come 
up with a new (!) nonlinear equation already on the second step. This equation has a 
non-trivial Backlund transform, admitting an interesting 27r-kink-shelf solution. 

The same technique was shown to be applicable to the study of considerably more 
difficult (l-|-2)-dimensional nonlinear PDEs. For instance, in [14] the proliferation proce- 
dure was successfully adapted to the Kadomtsev-Petviashvili and Boiti-Leon-Pempinelli 
equations. 

Let us apply this formalism to the HD equation. First note that the LA-pair for (19) 
can be written as a system of two Ricatti equations: 

9y = -Xg'^-Ug^l, 

(21) 

gt = x{2Uy-^-^-Ax) ^2-(^ + 4At/)^ + ^ + 4A. 

The second summand in the right hand side of the second equation has a term, denoted 
as Viy, representing a fairly long expression which can be derived from (19). The function 
g — g{y, t) is connected with the solution ip of (18) as ^ = ipy/ip- Excluding the function 
V from (1.1) and returning to the old variables via x = g, u = gy, we obtain a modified 
Harry Dym (mHD) equation: 

Ut = u^U3x + Su\xUxx - 3A\«2 ^ + — ^ + ^ 22 

X X'^ x-^ 
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(By analogy with the equation mKdV in [12], we call (1.2) the mHD equation.) As 

one can sec, this equation has a different countenance than the HD equation. However, 
omission of all the summands but the first one in the right hand side of (22), yields 
the HD equation (16) with f3 = 0. Note that (22) can be rewritten quite nicely in new 
variables x = 1/z, u{x, t) = \/0{z, t): 

{e-'/') ^ = ^-^ + z'(^lz%, + lz'd,, + 9z9, + 3d-3y (23) 

The formula (23) can be simplified even further by changing 

e = e-''^ri{^,t) + l, y = logz. 

As the result, it becomes 

However, we will be considering the mHD equation in the form (23). Note than in the 
stationary 9t = case, it reduces to a linear ODE! 

The dressing chain method produces not only the equation (23), but also its LA-pair. 
It is constructed as follows. Return to the chain (15) and let = 1/g, fn+i = ^, A„ = A, 
An+i = /U. Considering |U as a spectral parameter, one can see that (15) can be viewed 
as an L-equation of the LA-pair for equation (23). One should also define the second, 
non-stationary chain Ct for the functions gn (in terms of the dynamical equation in (18)) 
and build the A-equation. Omitting the lengthy but straightforward computation, we 
present the LA-pair for equation (23), written in the variables t, z: 



m = ^2 I 1 I _J, IlS] ^ _ 



(24) 



where 

a = 



z^ 



*t = +b'^ + c, 

4/x + 2A - 1^ - 2X^/e +{d-l)z^ + 2z^e, + \z^e,, 
b = 4:(^^-z-zV9^fx+lz^e,z + 2z'^9,+ {9-^9z,-l) 

3Xz^e, + 3X{i-e)z-^ + {^)', 

c = 4:ii-lz'^e,,-3z^e,+ (l-3e-2V9)z^-2X+{^)'^ . 

Note that the spectral parameter in (24) is /x, whereas A enters the non-linear equation 
(23). As one can see, the LA-pair for (23) also has the form of a pair of Ricatti equa- 
tions. These equations can be simultaneously linearized in order to represent (23) as a 
compatibility condition for two linear equations, as it is done in the theory of solitons. 

We will not proceed further with the mHD equation. To conclude this section, we 
would like to repeat the statement that the exact solutions of (23) are easily found via 
the dressing technique, and this procedure can be extended in order to produce the 
mHD equation and its LA-pair. It is worth emphasizing that the equations HD and (23) 
are members of different hierarchies. Thus, discrete symmetries enable one to establish 
connections between different integrable equation hierarchies, promoting the unification 
of knowledge about them. 
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MOUTARD TRANSFORMATIONS 



We devote this last section to equation (5) which has been obtained from the Maxwell 
equations (2) in the case of an isotropic but inhomogeneous in two directions (x, y) 
medium. 

Clearly, a PDE (5) is harder to investigate than an ODE (1). Nevertheless, its analysis 
in terms of the Darboux transform (10) known also as the Moutard transformation [6] 
is quite similar to its ODE cousin. Below we shall present the relevant formulae without 
the derivation details. 

Let i/j = iIj{x, y) and (j) = 4){x, y) be two particular solutions of (5), namely 

AV' - Xetl) ^ A(f)- Xe(f) = 0. (25) 

We choose the function ^ as a prop solution. Then the following transformations repre- 
sent an analogue of (10): 

V'^V^^) = %^, 6^6W=e-2AAln<^, (26) 


where 

e[ip, 4>]^ dx^e^, {<pd,ip - ^d,<j>) . (27) 

Above, the following (standard) tensor notations have been used: G {1, 2}, x^ G {x, y}, 
e^i, is a fully antisymmetric tensor with £12 — 1, summation is implied over 
repeated indices. Note that a one-form, which is being integrated in the formula (27) is 
closed in the case when and are solutions of (26). Hence, the shape of the contour 
of integration T in (27) is irrelevant. 

One can verify by direct substitution of (26,27) into the formulae (28) below, that the 
dressed function i/;^^) satisfies the dressed equation (25) (with the potential e^^\x^ y) and 
the same spectral parameter value A). 

The Moutard transformation (26) can be iterated several times, and the result can be 
expressed via Pfafhan forms [6] . Instead, we direct our interest to the Maxwell equations 

2 

(2). A straightforward computation (recall, X = —■^) yields the expressions for the 
dressed electric and magnetic fields E^-*^), B^^^: 



E(i) = e^-t (0,0, V^'^) , B(i) = ^e^-* V^^^), o) , = e(i)EW. (28) 

On the basis of (28) , one can build a variety of exact solutions of the Maxwell equations. 
As a simple example let's dress e = 0. This isn't quite a medium, but one can easily 
proceed with formal calculations (26-28) which result in a new "medium" whose dielectric 
permitivity e^^\x,y) and the stationary component of the field if)^^^ are as follows: 

^(1^ ^ 8c^ a'{z)h'{z) , (1) ^ a{z)(i{z) - a{z)h{z) + ^{z, z) 

^'^ {a{z) + h{z)f' a{z) + h{z) ' ' 



where 



^{z, z) = J dz {a{z)a'{z) - a{z)a'{z)) + J dz iP'iz)b{z) - b'{z)p{z)) , 
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a{z), a{z), b{z), (3{z) are arbitrary functions of z = x + iy, z = x — iy. Note that the 
function from (26,28) provides in fact a general solution of the dressed equation, 
for it is described in terms of two arbitrary functions a{z) and P{z). To ensure that the 
quantities found correspond to a physical non-absorbing medium, one should require that 
the dressed dielectric permitivity function e^^^ be real. This imposes an extra restriction 
to the quantities a{z) and b{z), namely b{z) = a{z). Generally speaking, the functions 
e*-^-* and i/j^'^^ from (29) will have singularities along certain curves in the (a;, y)-plane. 

The reflexionless B-potentials for the one-dimensional problem (1) above, possess point 
singularities on the real line (corresponding to zeroes of the function u{x)). Clearly, their 
2D-analogues, such as (29) for the equation (5) allow a much more diverse structure of 
singularities on the real plane. On the other hand, not requiring that the quantity e^-*^) be 
real, one obtains an absorbing medium which may not be devoid of interest for physical 
applications. 

In conclusion, let us study a dressing chain generated by the Moutard transforma- 
tions (26) . A simple periodic closing of the dressing chain results in a regular dielectric 
permitivity, similar to the ID case studied above. 

Denote /„ = ln0, fn+i = In'^*^-'^). Then after a straightforward computation 

A (/n + fn+l) = ||V/nf - ||V/n+if , (30) 

where || • || is the Euclidean norm. 

The chain (30) is closely related to that of Veselov and Shabat [10] for the Schrodinger 
equation. Choosing specifically as 

fn = \/Ky + J dxgn{x), 

and substituting it into (30) (A^ being constant), we obtain for the quantities gn{x) the 
following expression 

{gn + 9n+l) = 9n~ 9n+l + '^n " 

matching the corresponding formula of [10]. 

The simplest periodic closing of the dressing chain (30) is fn+i = fn = F{x, y), which 
implies that the latter function F is harmonic, and that the regular dielectric permitivity 
function in the corresponding medium is given by the formula 
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